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ABSTRACT 



Aims. The analysis of near infrared spectropolarimetric data at the internetwork at different regions on the solar surface could offer 
constraints to reject current modeling of these quiet areas. 

Methods. We present spectro-polarimetric observations of very quiet regions for different values of the heliocentric angle for the 
Fei lines at 1.56 yum, from disc centre to positions close to the limb. The spatial resolution of the data is 0.7 - 1". We analyze 
direct observable properties of the Stokes profiles as the amplitude of circular and linear polarization as well as the total degree of 
polarization. Also the area and amplitude asymmetries are studied. 

Results. We do not find any significant variation of the properties of the polarimetric signals with the heliocentric angle. This means 
that the magnetism of the solar internetwork remains the same regardless of the position on the solar disc. This observational fact 
discards the possibility of modeling the internetwork as a Network-like scenario. The magnetic elements of internetwork areas seem 
to be isotropically distributed when observed at our spatial resolution. 

Key words. Sun: magnetic fields — Sun: atmosphere — Polarization — Methods: observational 



1. Introduction 

The presence of magnetic fields i n the solar internetwork was 
discovered mor e than 30 years ago (Livings ton & HarvevLll975l 
Smit hsonlll975l) . Since that work, improvements on the instru- 
mentation have allowed the observation of the full Stokes vec- 
tor in those regions with a signal-to-noise ratio good enough to 
. £h ! retrieve the magnetic field vector from the observational data. 
. Nevertheless, the internetwork magnetic topology remains yet 
' very indistinct and sho ws a growing complexity as we improve 
j the quality of the data. Martini d 1987b observed that in her video- 
magnetographs at a spatial resolution of 3" the longitudinal 
component of internetwork magnetic fields was present every- 
where in the solar disc. She immediately concluded that these 
magnetic structures should have very small scales and a very 
tangled geometry, giving as an example a scenario in which the 
magnetic fields in the internetwork would consist of a maze of 
small loops. 

As spectro-polarimeters have made possible the precise 
detection of these signals, the efforts have concentrated on 
the study of the distribution of magnetic fields at disc cen- 
tre. Different methodologies (mainly based upon the use of 
the Zeeman and Hanle effects) shed light on different as- 
pects of the internetwork magnetism. Concerning the Zeeman 
effect, the Stokes Q, U and V signals detected on the 
most widely used spectral lines (Fei at 1.5 //m and 630 
nm) are very weak at the best spatial resolutions of 0.5- 
1". Moreover, the Stokes / profile seems to come from a 



fiel d free atmosphere, as expected from weak ly polarized me- 
dia ([Sanchez Almeida & Truiillo Buenol Il999h . The filling fac- 
tor of the magnetic elements retriev ed from the analysis of 
these signals is alway s around 2 % (Khomenko et al., 2003; 
ISanchez Almeida et all 120031; iMartmez Gonzalez et aUl2007al) . 
The rest of the resolution element would be filled with very 
weak magnetic fields. Another possibility is that the real ele- 
ment is filled with mixed polarity magnetic fields which would 
partially cancel out due to the lack of spatial resolution. The 
magnetic field strength distributions recovered for this 2 % of 
the resolution element appear to show a preference for magnetic 
fields around t he equipartition field at photospheric heights and 
even weaker (Khomenko et al., 2003; Martinez Gonzalez et al., 
| 2007at iRamfrez Velez & Lopez AristeL 120071) . Additionally, 
Asensio Ramos et al. (120071) have shown the first direct observa- 
tional evidence of flux cancellation in the internetwork (note that 
their internetwork region is surrounded by a very enhanced net- 
work structure), showing that more than 95 % of the magnetic 
flux is cancelled in the resolution element (~ 1"). This means 
that the above-mentioned distributions could not give a complete 
vision of the internetwork magnetism. lAsensio Ramos et all 
d2007l) give an amount of 250 G for the mean magnetic field 
in the resolution element (note that the values of the magnetic 
flux density obtained by means of the Zeeman effect are below 
10 Mx/cm 2 ). This would mean that the magnetism of the inter- 
network could play an important role on the solar global mag- 
netism . This has also been pointed out by works using the Hanle 
effect dTruiillo Bueno et aill2004l) . 



Send offprint requests to: M. J. Martinez Gonzalez 



2 



M. J. Martinez Gonzalez et al.: Near-IR internetwork spectro-polarimetry at different heliocentric angles 



The study of regions in different positions of the solar disc 
could represent a strong constraint to reject models for the so- 
lar internetwor k. Using the 630 nm lines with a spatial reso- 
lution of ~ 1", iLitesi Q2P02) built a histogram of both circular 
and linear polarization signals in two quiet regions, one at disc 
centre and another at // = 0.82 (being /i the cosine of the helio- 
centric angle). He did not observe significant linear polari zation 
signal s in neither of the two regions. However, Figure 9 in lLitesI 
(2002) shows that the histograms of circular polarization do not 
present any va riation for those signa ls whose integrated signal 
is below 0. 005.lMeunier et all ( 19981) . studying integrated polar- 
ization, and H arvey et al. (2007) using magnetograms, show the 
presence of a horizontal component of the magnetic field every- 
where in the solar disc. In this paper we present the first study of 
the internetwork at several positions on the solar disc using high 
quality 0.8" spectro-polarimetric data. 



2. Observations and data reduction 

The observations consist of 2-dimensional maps of the very quiet 
Sun taken at different positions covering the two solar hemi- 
spheres. In each pixel of these maps we recorded the four Stokes 
parameters of the Fei lines at 1.56 /mi. The observations were 
performed at the Vacuum Tower Telescope (VTT, Observatorio 
del Teide) during August 2000 and Jul y 2006 with the Tenerife 
Infrared Polarimeter (TIP) instrument (Coll adoslll999l) . The ob- 
servations during 2006 made use of the adaptive optics system 
attached to the telescope. This led to a much improved spatial 
resolution and to an important decrease in the exposure time 
required to reach a noise level comparable to those obtained 
for the data of August 2000. The maps of August 2000 cover 
the following positions (we also indicate the quadrant of each 
scan): //=0.88 (N), //=0.4 (N) and /i=0.28 (E). The maps of July 
2006 were taken at ju=l (disc centre), //=0.81 (SW), //=0.73 (S), 
yi/=0.62 (W) and //=0.46 (W). Active regions or very enhanced 
Network were avoided using the Ca II K filter available simul- 
taneously to the observations, since bright regions in this filter 
are go od indicators of magnetic flux concentrations dLites et al.L 
119991) . 

The data reduction consisted in the subtraction of dark 
current, flatfield correction and demodulation of the images. 
After that, we also corrected for o ther residual patterns, 
mostly polarization-dependent fringes dSemell [2003). Most of 
the instrumental crosstalk can be removed from the data us- 
ing the calibration optics locate d before the beam splitter 
dSchlichenmaier & Colladosll2002h . However, the coelostat con - 
figuration of the telescope has to be modeled dColladosl fl999). 
The residual crosstalk from Stokes / to Stokes Q, U and V was 
removed by forcing the continuum of the polarization profiles 
to zero. The residual crosstalk between Stokes Q, U and V is 
very difficult to remove and a few percent may still remain. The 
last step of the data reduction was a de-noising procedure based 
on Principal Component Analysis. The noise level in the polar- 
ization profiles of the 2006 data sets is in the range 3 x 10~ 5 - 
8 x 10~ 5 in units of the mean continuum intensity, I c . The value 
in the 2000 data is ~ 10~ 4 I c . This means that the signal to noise 
ratio is ~ 15 (taking the most probable polarization amplitude as 
8 x 10~ 4 Ic) and ~ 8 in the 2006 and 2000 data sets, respectively. 
The spatial resolution is ~ 1" for the 2000 data and 0.7 - 0.8" 
for the 2006 observations. 

In order to analyze the data in a reliable way, we select those 
profiles which have a degree of polarization Ap (maximum of 
Vfi 2 + u 2 + v 2 ) higher than 4 x 10 4 I c . Following this crite- 



rion we do not introduce much bias in the results since we se- 
lect all kind of magnetic field vector configurations. In the disc 
centre observations we have detected some small accumulations 
of high Stokes V amplitude (> 0.01 I c ) that may correspond to 
Network points. Far from the disc centre one expects Network 
fields to be more inclined. Then, we reject those points with am- 
plitudes of total polarization higher than > 0.01 I c . This criteria 
is not perfect since linear polarization does not grow as fast as 
circular polarization decreases when the magnetic field inclina- 
tion with respect to the line-of-sight increases. Anyway, these 
points represent at most 1 % of the field of view. The selected 
internetwork profiles account for at least 95 % of the observed 
field of view in all maps. 



3. Analysis of circular and linear polarization 

From the selected profiles we compute the amplitude of the cir- 
cular polarization (Ay) in the following way. We find the position 
of the maximum and the minimum of the Stokes V profile. We 
compute the mean value of the profile around these positions (we 
take 3 points before and after the maximals). The amplitude Ay 
is the mean value of the absolute value of the amplitude of the 
two lobes. The linear polarization (Ai) is defined as the mean 
value around the maximum of Vfi 2 + U 2 . The same calculation 
is performed to obtain the total polarization Ap. Both the 1.5648 
and 1 .5652 /mi spectral lines show the same behaviour and we 
present here results for just the 1.5648 /mi line whose larger am- 
plitudes in the Stokes parameters makes the results less subject 
to noise effects. 

The comparison of the histograms of the polarization ampli- 
tudes between the data of 2000 and that of 2006 proved to be 
very difficult. First, the different spatial and spectral resolutions 
(the sampling in wavelength in 2000 data was twice the one of 
the 2006 data) and different noise levels result in a different po- 
larimetric sensitivity. Second, both observations have different 
number of pixels. This can be partially solved by normalizing 
all the histograms to its area. However, the fact that both data 
sets have different polarimetric sensitivities makes its compari- 
son impossible after the normalization. Figure Q] shows the de- 
gree of polarization {Ap) of the two different data sets. There 
is no evident difference between the signals recorded either on 
2000 or on 2006. However, for the sake of having reliable con- 
clusions we only explicitly compare the histograms for the 2006 
data. In order to compare all the observations we analyze other 
properties of the histograms that do not depend on the normal- 
ization. 

The top and center panels of Fig. [2] show the histograms of 
the circular and linear polarization amplitudes of the 1 .5648 /mi 
line for all the observed positions on 2006. We have normalized 
the Ay and Ai values to the mean continuum intensity of each 
map to avoid the effect of limb darkening. There is no evident 
variation of both circular or linear polarization signals in all the 
studied positions. However, the histogram of linear polarization 
at /(=0.62 is different from the others. A visual inspection of the 
data reveals that this data set is highly contaminated by cross- 
talk from Stokes V to Q and U. We conclude that, irrespective 
of the position on the solar disc, the polarimetric signals on the 
internetwork seem to be the same at our 0.7 - 0.8" spatial reso- 
lution. The bottom panel of Fig. [2] shows that the ratio between 
the circular and linear polarization does not change significantly 
over the solar disc. Since the signal to noise ratio is at least 15, 
we do not expect the results to be dominated by noise effects so 
that we can affirm that the observed pattern is of solar origin. 
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Fig. 1. Histograms of the total polarization Ap for all the ob- 
served areas on August 2000 (those corresponding to /i = 
0.88,0.4,0.28) and July 2006 (those corresponding to fi = 
0.8 1 , 0.73, 0.62, 0.46). The improvement of the spatial resolution 
from 1 to 0.7-0.8" has made it possible to detect a large amount 
of very weak signals. The results show that there is no varia- 
tion of the signal for different heliocentric angles. The vertical 
dashed line indicates the chosen detection limit. 

According to the bottom panel of Fig. [2] the most prob- 
able value for the ratio between circular and linear polariza- 
tion is close to 1 . This can be a surprising result if we have in 
mind all the results found in the literature concer ning the inter - 
network magnetism using the 630 nm lines (e.g., iLitesl 12002b . 
At this spectral range n o linear polarization can be found at 
0.5-1" spatial resolution dLitesl l2002t iMartfnez Gonzalez et al.L 
2007b). However, the linear polarization sensitivity of the 1.5 
pm pair of lines is very different from th at of the 630 nm 
pair. [Landi de glTnnocenti & Landolfil d2004l) define some use- 
ful magnitudes that can be used to get an idea of the sensitivity 
of a spectral line to the circular and linear polarization. They are 
defined as: 



s v = 



l ref 



sl = (- — ) Gd c , 

Arej 



(1) 



where s-y is the circular polarization sensitivity index and Si 
is the corresponding linear one, while A re f is a fixed reference 
wavelength. Consequently, these numbers have only sense when 
comparing two spectral lines. The symbol g stands for the effec- 
tive Lande factor of the transition and d c is the central depression 
in terms of the continuum intensity, defined as: 



d c - 



(2) 



being I(Aq) the intensity at the line core. The symbol G plays the 
same role as the effective Lande factor but for the linear polar- 
ization and it is defined as: 



(3) 



where 6 is a qua ntity that depends on the quantum numb ers of 
the transition (see lLandi deglTnnocenti & Landolfill2004h . Both 
1.5648 fi and 630.2 nm spectral lines have 5 = 0. The ratios be- 
tween the circular and linear polarization indices of the infrared 
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Fig. 2. Histograms of circular polarization (upper panel), linear 
polarization (central panel) and its ratio (bottom panel) for all 
the observed areas on July 2006. There is no variation of the 
distributions with heliocentric angle. 



and visible spectral lines are: 

(Sv)l.5imt 



(Sv)630nm 
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1.37 



= 4.17. 



(4) 
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This means that the circular polarization sensitivity is similar in 
both spectral ranges but the infrared lines present a linear polar- 
ization sensitivity that is 4 times larger than the visible lines. This 
fact explains the observed behaviour of the presence of ubiqui- 
tous linear polarization signals in the infrared and the lack of 
them in the visible for similar spatial resolutions. 

All the histograms of Figs. [TJ and [2] (and also the circular and 
linear polarization histograms of year 2000) present a common 
behaviour: a peak is found for values markedly above the detec- 
tion limit. In the case of the total polarization, we can see that the 
maximum of all histograms is far from the chosen noise thresh- 
old of 4 x 10~ 4 I c . Also the drop towards small amplitudes in 
the case of the circular and linear polarization occurs for values 
larger than three times their noise level. Such behaviour can be 
expected from cancellations of magnetic fields due to a decrease 
of the spatial resolution. In any case we do not discard a bias 
introduced by noise or other systemathic defFects (e. g. interfer- 
ence fringes) in the determination of the amplitudes. This issue 
is now under investigation. 

In order to further compare the whole data sample we choose 
two different functional forms to fit the tails of the histograms. 
First, a combination of two power laws (y oc xr'). The first 
one (parameterized by b\) describes the drop from the peak to 
~ 7 x 10~ 3 I c for circular polarization or to ~ 3 x 10~ 3 I c for 
linear polarization, while the second one (parameterized by bj) 
describes the drop for the rest of the tail above these points. 
Second, a decreasing exponential law (y oc exp[-jt/a]) start- 
ing from the peak of the histograms. This second option is cho- 
sen following previous works on quiet Sun magnetism in the 
infrared. These parameters (a and bi) are independent of the nor- 
malization of the histograms and they are useful to compare the 
behaviour of these histograms with the heliocentric angle. 

Figure [3] shows the dependence of these three parameters 
with the heliocentric angle. Error bars represent statistical un- 
certainty. The left panels present the results for the circular po- 
larization amplitudes while the right panels refer to the linear po- 
larization amplitudes. No clear trend is found. Perhaps, it could 
be possible to identify a trend in the value of bt with fi, indi- 
cating the presence of magnetic fields with a preferential ori- 
entation. In any case, the poor statistics on the tails make this 
conclusion delicate. We cannot be sure if this behaviour is due 
to some isolated Network patches, high magnetic flux pixels that 
may contaminate the histograms or just statistical noise. In any 
case they represent a small percentage of the points under study. 
The key fact is that the majority of the signals (coming from in- 
ternetwork areas) have no trend with the heliocentric angle. This 
analysis states that all the positions on the internetwork quiet Sun 
are indistinguishable from the point of view of the distribution 
of polarization signals. Our data at 0.8 - 1" reveals no clear vari- 
ation of the internetwork magnetism with the heliocentric angle, 
supporting the idea of an isotropically distributed field. 

4. Analysis of amplitude asymmetries of the Stokes 
V profiles 

For the computation of the Stokes V asymmetries we consider 
only the regular (two-lobbed) profiles. In order to select them, 
we first determine the number of lobes (and their position in 
wavelength and amplitudes) of both the 1 .5648 and 1 .5652 yum 
lines. Due to small residual patterns of the reduction, we can 
find some profiles where one of the spectral lines presents two 
lobes and the other one presents three. In order to extract the real 
lobes (whose number should be the same in both lines) we im- 
pose both spectral lines to have zero-crossings compatible with 



an error of 0.6 km/s. If the line has two lobes, the zero-crossing 
is defined as the point where the Stokes V line profile reaches 
zero in between the lobes. If it presents three lobes there will be 
two zero crossings: one between the red and the centre lobes and 
a second one between the blue and centre ones. Finally, we only 
select the two lobes that are compatible. This data selection re- 
duces the sample of "regular" profiles to 63+3 % of the field of 
view, which still represents a large fraction and the results will 
have a strong statistical significance. 

Figure [4] shows the am plitude asymmetry of th e regular 
Stokes V profiles, defined as Solanki & Stenfl 



try oi tn e 
ol d 1984 : 



CLb — Cl r 

ab + a r 



(5) 



where 6a accounts for the amplitude asymmetry while aj and 
a r refer to the absolute values of the blue and red lobes of the 
Stokes V profile, respectively. The area asymmetry is computed 



as: 



SA 



A r 



A b +A r 



(6) 



being SA the area asymmetry while Ab and A r are the abso- 
lute values of the areas of the blue and red lobes, respectively. 
Computing the area asymmetry implies choosing the boundaries 
for the integration of each lobe. In our case the integration in- 
terval is different for each profile. We select the initial point for 
the integration of the blue lobe as the position on the blue wing 
closest to the zero crossing which has a value lower than 3xl0~ 5 
I c while the ending point is chosen to be the zero crossing of the 
profile. This last one is also the lower limit for the integration 
of the red lobe and the final one is also the closest point to the 
profile with a value lower than 3xl0 -5 I c but in the red wing. 

As shown in Fig. [4] the amplitude and area asymmetries do 
not change significantly in all the studied positions. Small vari- 
ations can be detected but no clear trend with fi is found. The 
results for t he amplitude asymmetr y are compatible with those 
obtained bv lKhomenko et all d2003l) at disc centre, including the 
upper limit of 0.6 for the positive values. For the area asymmetry 
we found a wider range of values even though the distributions 
are symme tric and centered at SA = 0, also consistent with what 
is found by Khomen ko et al.l (|2003). 

Figure|5]shows the Stokes V amplitude versus the amplitude 
(left panel) and area asymmetries (right panel). Left panel shows 
that negative amplitude asymmetries appear to present ampli- 
tudes preferentially in the range between 3 x 10~ 4 and 4 x 10 -3 
I c . For the positive asymmetries the case is more complicated, 
where one is able to distinguish at least three different regimes. 
Like for the case of negative asymmetries, Stokes V amplitudes 
in the range 0.3 - 4 x 10~ 3 I c seem to cover the whole range of 
values (up to an upper limit of 0.6). For higher amplitudes, the 
higher the value of Ay the smaller the value of 5a. Interestingly, 
amplitudes lower than ~ 3 x 10 -4 I c show the inverse case: the 
higher the amplitude the higher the asymmetry. 

5. Discussion and conclusions 

We have analyzed high quality spectro-polarimetric data of the 
Fe i lines at 1 .56 fim in order to infer information about the in- 
ternetwork magnetism at different positions on the Sun's surface. 
The whole field of view presents significant signal, meaning that 
the magnetic fields pervade the observed areas. 

We have found that the circular and linear polarization am- 
plitudes do not have any clear dependence on the heliocentric 
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Fig. 3. Upper panels: exponents of the two power laws {b\ shown in stars and bi shown in squares) that fit the tails of the distributions 
of circular and linear amplitudes shown in Fig. [2] Bottom panels: values of scale factors a for the tails of the distributions. 



K3 



10.000 I 



1.000 



0.100 = 



o 



0.010 z 



0.001 




10.00 4 



1.00 



0.10 



0.01 



0.6 




Fig. 4. Histograms of the amplitude (left panel) and area (right panel) asymmetries for the regular Stokes V profiles. 



angle. This fact goes against a Network-like scenario for the 
internetwork: quasi-vertical flux tubes cannot explain this ob- 
servational result, nor in fact any field topology with a pre- 
ferred orientation within the field-of-view. An isotropical dis- 
tribution of magnetic fields, oriented in all directions in the 
whole field of view, is on the other h a nd exp ected to show 
this behaviour. iMartfnez Gonzalez et al.l (|2007c) found that at 



least 10-20 % of the magnetic flux in the internetwork is con- 
ne cted by l ow-lying loops. Consequently, the scenario proposed 
bv iMartinl (1 19871) of an internetwork characterized by a myr- 
iad of small loops is a very reasonably idea that is compatible 
with all the observational constraints presented in this work. 
Of course one can think of other scenarios that are compat- 
ible with the observations. (Stenfl^, 119871: iManso Sainz et all 
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Fig. 5. Amplitude (left panel) and area (right panel) asymmetries plotted versus the amplitude of the regular Stokes V profiles. 



I2004t iTruiillo Bueno et al.L [2004 adopt tur bulent internetwork 
magnetism. ISanchez Almeida & Litei (2000) proposes a micro- 
structuration of the atmosphere (MISMA) to explain all the mag- 
netic phenomena on the solar surface. All of them are compat- 
ible with the presented results and our efforts should be headed 
towards finding more constraints to reject some of them and 
strengthen others. 

The size of the magnetic structures can also be constrained 
by the information presented in this study. First, improving the 
spatial resolution, we do not s ee a global increase in the signals. 
iLites & Socas-Navarrol (|2004|) found no incr ement of the mag- 
netic flux density from 1" to 0.6". Recently ILites et all (2007) 
computed a magnetic flux density of about 1 1 Mx/cm 2 at 0.3" 
using HINODE's d ata, wh i ch is c ompatible with the value of 10 
Mx/cm 2 found bv lMartinl d 19871) at 3". This means that, either 
the magnetic field structures are already resolved at ~ 0.5" or we 
are very far from resolving them. The fact that the polarimetric 
signals do not vary along the solar surface would point towards 
very small structures as the responsibles for the internetwork 
magnetism. The size of these magnetic structures is something 
not yet constrained by the observations. 

We have presented a study of high quality spectro- 
polarimetric data in different positions of the solar surface, from 
the disc centre towards ju = 0.28. This is the first step of the 
study of the variation on the magnetism of the internetwork with 
the heliocentric angle. Much more work has to be done by re- 
trieving physical parameters as the magnetic field strength vec- 
tor, magnetic flux, etc. to really constrain the modeling of the 
internetwork and reject models that are not compatible with the 
results. 
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